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I The majority of nearby early-type galaxies contains detectable amounts of emission-line gas 

■<!;;j- ■ at their centers. The emission-line ratios and gas kinematics potentially form a valuable 

I diagnostic of the nuclear activity and gravitational potential well. The observed central gas 

■ velocity dispersion often exceeds the stellar velocity dispersion. This could be due to either 
the gravitational potential of a black hole or turbulent shocks in the gas. Here we try to 

_ discriminate between these two scenarios. 

O ' 1.1 Introduction 

> ' Optical emission-line gas is detected in at least 50% of the centers of nearby early- 

^ [ type (E/SO) galaxies (e.g., Phillips et al. 1986; Goudfrooij et al. 1994; Macchettoet al. 1996; 

Ho et al. 1997). The gas velocity dispersion dgas (corrected for instrumental broadening) 
typically increases towards the nucleus, reaching values similar to the the central cTstar as 
' measured through 100s of pc scale sized apertures (e.g., Zeilinger et al., 1996). Closer to the 

^ . nucleus, Cgas often exceeds tTstar (e-g-, Ferrarese et al. 1996; Macchetto et al., 1997; van der 

Marel & van den Bosch 1998; Cappellari et al. 2002; Verdoes Kleijn et al. 2002). Potential 
contributors to o-gas are motions due to (i) a black hole (BH) gravitational potential or (ii) 
shocks/turbulence in the gas. In the former case, dgas could be caused by three dimensional 
motion of collisionless cloudlets of gas. More likely, the collisional gas is expected to have 
settled in a disk. In that case, cTgas will be caused by differential rotation over the aperture. 
However, shocks and/or turbulence can contribute to Ugas if there is a permanent input of 
kinetic energy (perhaps the active galactic nucleus) to sustain them. Thus, determining 
the driver of cTgas has implications for measurements of the BH mass (Mbh) and the flow 
properties of potential BH accretion material. Here we present preliminary modeling of dgas 
and (Tstai in the nuclei of 16 active and 4 quiescent early-type galaxies to constrain the relative 
importance of gravitation and shocks/turbulence. 



1.2 Sample and Data 

The galaxy sample consists of 16 active galaxies with radio-jets taken from the 
UGC FR I sample (see Noel-Storr et al., these proceedings, for selection criteria) and four 
relatively quiescent early-type galaxies (NGC 3078, NGC 4526, NGC 6861 and IC 1459) 
(Table 1). The four quiescent early-type galaxies were selected to have central dust and/or 
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Table 1.1. Early-Type Galaxy Sample. 



Galaxy 


Type 


Mb 
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Type 


Mb 


D 








(Mpc) 








(Mpc) 


(1) 


(2) 


(3) 


(4) 


(1) 


(2) 


(3) 


(4) 


NGC315 


E+: 


-21.9 


67.9 


M87 


Eh-0-1 pec 


-21.4 


15.4 


NGC 383 


SAO-: 


-21.4 


65.2 


NGC 5127 


E pec 


-20.6 


64.4 


NGC 541 


cD 


-21.2 


73.3 


NGC 5490 


E 


-21.1 


69.2 


NGC 741 


EO: 


-21.9 


70.2 


NGC 7052 


E 


-19.9 


55.4 


3C 66B 


E 


-21.1 


84.8 


3C 449 


SO-: 


-19.6 


68.3 


NGC 2329 


SO-: 


-21.1 


76.7 


NGC 7626 


E pec: 


-21.1 


46.6 


3C 264 


E 


-21.0 


84.4 


NGC 3078 


E2-3 


-20.6 


35.2 


UGC7115 


E 


-20.4 


90.5 


NGC 4526 


SAB(s)0 


-20.5 


16.9 


NGC 4335 


E 


-20.3 


61.5 


NGC 6861 


SA(s)0-: 


-20.2 


28.1 


M84 


El 


-20.8 


15.4 


IC 1459 


E3 


-21.1 


29.2 



The galaxies are taken from the UGC FR I radio galaxy sample (Noel-Storr et al., 
these proceedings) except for NGC 3078, NGC4526, NGC 6861 and IC 1459. 
Col. (2): Hubble classification from NED. Col. (3): Absolute blue magnitude from 
LEDA. Col.(4): Distances from Faber et al. (1989), Tonry et al. (2001), or, if not 
available, directly from recession velocity and Hq = 75kms"'Mpc"'. Col. (5): 
Galaxy is part of the UGC FR I sample (1) or not (0). 

gas disks. Four of the 21 UGC FR I galaxies were not included here because either dstar is 
not available or the location of the nucleus could not be determined reliably, mainly due to 
dust obscuration 

HST/STIS (HST/FOS for IC 1459) emission-line spectra are available for the sample. 
They include the HaH-[NII] and [SII]6716,6731 lines. A central dgas per galaxy is deter- 
mined by fitting single Gaussians to each of the lines in Hq!H-[NII] and [SII]6716,6731 at the 
nucleus. An analytic emission-line flux profile is determined by fitting a double exponen- 
tial to the emission-line surface brightness measurements from the HST/STIS spectra and 
HST/WFPC2 emission-line images (where available). The central stellar dispersion cTstai for 
the sample galaxies was obtained from the LEDA catalog*. 

1.3 CoUisionless Gas Models 

Can the gravitational potential of a BH mass induce gas motions which explain 
fjgas? Current observations are consistent (e.g., Kormendy & Gebhardt 2001, for a review) 
with most (and perhaps all) early-type galaxies containing a black hole with a mass related 
to (Tstai, i.e., the 'MBH-cTstar' relation (Gebhardt et al. 2000; Ferrarese & Merritt 2000). We 
compute the expected Ugas for such a BH mass under idealized circumstances. We assume 
the gas is a collisionless isotropic spherical distribution of gas cloudlets orbiting the BH. 
The predicted gas velocity dispersion o'gas(sphere,l) takes into account aperture size, PSF 
and instrumental line width. 

Fig. ll.lh shows that Ugas typically exceeds cTgas (sphere, 1). However, the calculation ne- 

* LEDA database can be found at |http : / / leda ■ univ-lyonl . f r| 
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Fig. 1.1. The observed gas velocity dispersion (Xgas as a function of predicted ve- 
locity dispersions for radio galaxies (open circles) and relatively quiescent galaxies 
(filled circles), (a): The model for (Tga,s(sphere,l) assumes an isotropic spherical 
distribution of gas cloudlets orbiting a BH with a mass according to the Mbh - cr^xm 
relation, (b): similar to (a) but now including the stellar mass distribution as well. 
The dashed and solid lines indicate o'gas(sphere)= (l-y/1/2, 1,^/2) x a^^^ respec- 
tively (to facilitate comparison to the Mbh - CTstar)- The figure suggests the gas 
might be largely coUisionless: see text for details. 



gleets the stellar mass present, which will increase (Tgas(sphere,l). To illustrate the possible 
effect of this. Fig. II. lb shows crgas(sphere,2) assuming an isothermal spherical stellar mass 
distribution in addition to the BH mass, i.e., o'g25,(sphere,2) = crg^s(sphere, \ ) + a^^^ typ- 
ically agrees within a factor 2 to a^^J^sph&rt,!). In other words, the Mbh needed to account 
for cTgas agrees within a factor ^ 2 with Mbh as predicted by the Mbh - Cstai- relation (because 
Mbh ^ cTaas)' under these simplifying assumptions. The estimated intrinsic scatter in the 
relation is also a factor 2 (Tremaine et al. 2002). Thus, an additional contribution to o-gas 
from collisions (e.g., shocks/turbulence) is not required in this case. 

Alternatively, we can assume that the collisional gas settles in a disk. We idealize the disk 
as an infinitely flat disk of gas particles/cloudlets, inclined by 60° to the line of sight, and in 
circular rotation in the same combined BH and stellar potential as used in the spherical case. 
On the nucleus, the differential rotation will broaden the emission-line flux profile. Fig. ll.2h 
shows that the predicted gas velocity dispersion CTgas(disk) roughly agrees with a„.^^. Hence 
no shocks/turbulence are required on average; a similar result as for the gas spherical model. 

Gas dynamical models of galactic nuclei, based on gas mean velocity Vgas instead of o-gas 
often reach a different conclusion. For five sample galaxies, the central BH mass has been 
measured directly from gas mean velocities Vgas in extended gas disks. Fig. 11.2b shows 
that the velocity dispersion crgas(vgas) predicted by our gas disk model for these BH masses 
underpredicts Cgas in four cases. Hence, an additional contributor to the velocity dispersion 
is needed to explain the observed values in our gas disk model. This agrees with the mean 
velocity modeling which invokes an additional source of dispersion, 'gas turbulence', to 
explain dgas (Maciejewski & Binney 2001 ; van der Marel & van den Bosch 1998; Cappellari 
et al. 2002; Verdoes Kleijn et al. 2002). Only for NGC 4486 does a detailed gas disk 
modeling by Macchetto et al. (1997) account for both Vgas and Cgas, consistent with our 
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Fig. 1.2. The observed gas velocity dispersion (Xgas as a function of predicted ve- 
locity dispersions for radio galaxies (open circles) and relatively quiescent galaxies 
(filled circles), (a): The model assumes an infinitely thin gas disk in circular ro- 
tation, inclined by 60° to the line of sight and a BH mass Mbh according to the 
MBH-CTstar relation, (b): same as (a), but now for Mbh inferred from published, 
detailed modeling of the gas mean velocities Vg^j. The figure suggests that Vgas and 
(Tgas cannot be simultaneously fit by the gas disk model: see text for further details. 



modeling. It is reassuring that our simplistic modeling is consistent with the more detailed 
gas dynamical models on the difference between NGC 4486 and the other four galaxies. 

It is less reassuring that Fig. II. 2b suggests that a thin, circular gas disk model in general 
cannot explain simultaneously observed gas velocities and dispersions. For a simultaneous 
fit one has to invoke gas turbulence. Alternatively, the extended gas might somehow rotate 
at subcircular velocities, which would lead the models to infer a lower BH mass and hence 
lower nuclear differential rotation. 

Qualitatively, Fig. ll.ll suggests another possible solution to this problem: the gas changes 
from a flat disk configuration outside the nucleus to a (more) spherical distribution closer to 
the nucleus. The nuclear spherical gas distribution explains the central o-gas. Furthermore, 
the new configuration implies that (i) the flat gas disk effectively has an inner radius and (ii) 
only part of the observed nuclear emission-line flux is coming from a disk. Both changes put 
the gas associated with the circular disk at a larger mean distance from the BH, compared to 
the continous disk model, and hence at lower circular velocities for the same BH mass. 

1.4 Discussion & Conclusions 

The simple models discussed above suggest that the observed excess in velocity 
dispersion in nuclear emission-line gas disks could be explained by a more spherical gas 
distribution towards the nucleus. In such a model, the gas remains collisionless (as in the 
customary single gas disk model) and thus avoids invoking ad hoc turbulence or shocks 
which require a permanent input of kinetic energy to persist. Nevertheless, it could well 
be that activity associated with the BH can provide such a continuous input of energy. The 
outliers in Figs. ll.ll & ll.2l and the tendency for the collisionless models to underpredict Cgas 
more for active galaxies rather than for quiescent galaxies perhaps favor this scenario. 
More sophisticated modeling and further observations can help to determine between 
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different gas configurations. On the one hand, more accurate modeling of the stellar mass 
density and the emission-line velocity profiles are expected to decrease the predicted Cgas 
at fixed BH mass. The central density profile of bright ellipticals is generally shallower 
than an isothermal distribution. Emission-lines often show extended wings: the true 2nd 
moment is then larger than the Gaussian 2nd moment. This effect is taken into account in 
the gas disk model, but not in the spherical model. On the other hand, emission-Une flux 
profiles are typically marginally resolved. They may contain a significant contribution from 
an unresolved component which can be arbitrarily close to the BH. This might increase the 
predicted a^^^. 

UV spectroscopy can determine directly if shocks are present. In the presence of shocks 
the gas can be shock-ionized. If absent, the emission-line gas is only photo-ionized. The two 
ionization mechanisms produce different UV emission-line ratios (e.g., Dopita et al. 1997). 

In summary, the current results cannot distinguish between the various gas configurations. 
However, accurate BH mass measurements (and hence BH demography) and BH accretion 
models depend on such a distinction. Improvement of the simplistic gas dynamical modeling 
presented here and UV observations can help to resolve this. 
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